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Abstract 

We investigate the signatures of Fermi liquid formation in the M = 4 super Yang-Mills 
theory coupled to fundamental hypermultiplet at nonvanishing chemical potential for the 
global U(l) vector symmetry. At strong 't Hooft coupling the system can be analyzed in 
terms of the D7 brane dynamics in AdS$ x S 5 background. The phases with vanishing and 
finite charge density are separated at zero temperature by a quantum phase transition. 
In case of vanishing hypermultiplet mass, Karch, Son and Starinets discovered a gapless 
excitation whose speed equals the speed of sound. We find that this zero sound mode 
persists to all values of the hypermultiplet mass, and its speed vanishes at the point of 
phase transition. The value of critical exponent and the ratio of the velocities of zero and 
first sounds are consistent with the predictions of Landau Fermi liquid theory at strong 
coupling. 
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1. Introduction and summary 

Gauge/string duality is a spin-off of string theory which allows access to the dynamics 
of strongly coupled field theories by relating them to string theories in certain backgrounds. 
In principle, this relation provides a unique way of understanding phenomena at strong 
coupling which cannot be analyzed by other means. One of the fundamentally interesting 
questions (with many practical applications) concerns the behavior of the Fermi systems at 
strong coupling. Here we attempt to analyze the possibility of the formation of Fermi liquid 
in a very specific system - Nf Af = 2 fundamental hypermultiplets coupled to SU(N C ) 
M = 4 super Yang Mills theory. As usual, gauge/string duality methods are applicable in 
the regime of Af— >oo, Nf finite, and large four-dimensional 't Hooft coupling A. Recent 
work on this system in the strong coupling regime includes [1-13]. 

The hypermultiplet field content includes two Weyl fermions and two complex scalars. 
The global symmetry of the theory is U(Nf), but we will be mostly concerned with its 
U(1)b subgroup 1 and our results are not going to depend on the value of Nf. We will be 
interested in the situation where nonvanishing chemical potential for U(1)b is turned on. 
In the case of free fermions, such a setup would result in the formation of a degenerate 
gas of fermions, whose Fermi energy is equal to the value of the chemical potential ji. 
As the interaction is turned on, possible scenarios include formation of the Fermi liquid, 
superfluid, superconductor, or, possibly, some other exotic phase. 

Fermi liquid is a state whose properties can be described in terms of the dynamics of 
quasiparticles of Fermi statistics taking place in the narrow region around Fermi surface. 
The momentum and energy of a quasiparticle in the degenerate Fermi liquid at T = 
are bounded from above by the Fermi momentum kp and Fermi energy. Non-relativistic 
Fermi liquids are characterized by a number of properties. Among them are heat capacity 
proportional to T and the existence of zero sound. In [14] the D3-D7 system with vanishing 
hypermultiplet mass m has been investigated. It has been noted that the heat capacity of 
the system is proportional to T 6 (unlike the normal Fermi liquid). However, on the basis 
of the existence of a gapless excitation (arising as a massless pole in the density-density 
correlator) [14] argued that a novel kind of quantum liquid might be formed. This massless 
excitation is associated with zero sound in [14] , and we discuss some issues related to this 
identification in this paper. 

1 The subscript here stands for "baryon symmetry" . 
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We generalize the work of [14] to the finite value of to and find that the gapless 
excitation persists and the dispersion relation has an interesting dependence on to. In 
the relativistic regime, \x to, we reproduce the results of [14] where the speed of zero 
sound uq was found to be equal to the speed of regular sound. On the other hand, in 
the non-relativistic regime, ^ - m € m, we find that uq is proportional to {p, — to) 2, and 
differs from the value of the first sound, u\ by an overall multiplicative constant. Such a 
behavior is consistent with the prediction of Landau Fermi liquid theory. In particular, 
the value of the critical exponent is precisely the expected 1/2, while the ratio of u$ and 
u\ is expected to go to a constant in the limit of strong coupling. 

The rest of the paper is organized as follows. In the next section we give a brief 
review of the relevant parts of Landau's theory of Fermi liquids. In Section 3 we review 
the thermodynamics of the D3-D7 system. At large value of the 't Hooft coupling the 
dynamics of fundamental matter is captured by the DBI action for D7 branes propagating 
in the AdS§ x S 5 background [15]. At low temperatures, increasing the value of chemical 
potential results in a phase transition between the phases with vanishing and non- vanishing 
charge density. In the brane setup this corresponds to the brane embedding intersecting 
the horizon becoming thermodynamically preferred. In the limit of small temperature, the 
thermodynamics can be studied analytically [16]. 

In Section 4 we derive equations of motion for the fluctuating fields on the brane. We 
largely follow [17] and [15], with new ingredient being the finite value of the hypermultiplet 
mass to. We then proceed to analyze the equations in the regime of small frequency and 
momentum, find the massless excitation and compute its velocity as a function of the 
hypermultiplet mass to. We discuss our results in Section 5. Some technical results appear 
in the appendices. 

2. Landau Fermi liquid theory 

In this section we review the basics of Landau Fermi liquid theory and quote the results 
for the velocities of first and zero sound. We also show that the ratio of these velocities goes 
to a constant for generic strong interactions between the quasiparticles. Non-relativistic 
Fermi liquid at zero temperature involves the dynamics of quasiparticles of Fermi statistics 
whose dispersion relation is not necessarily equal to that of free fermions. The momentum 
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distribution is described by a step function 6>(|k| — kp) where kp is the value of the Fermi 
momentum, which is related to the particle density via 

N k% , rt „, 
n = — = — o" 2.1 
V 37t 3 h s 1 ' 

The quasiparticle description is, in fact, only valid in the small vicinity of the Fermi surface. 
The change in the total energy of the system due to the change in the distribution function 
5n(k) is 

5E = J e(k)5n(k)dr + J /(k, k')8n{k)5n{k')dTdT' (2.2) 

where dr = Vd 3 k/(27rh) 3 is the element of the phase space, e(k) is the energy of a quasi- 
particle with momentum k (whose absolute value we denote by k). In eq. (2.2) the second 
term describes interaction of the quasiparticles which takes place only in the narrow re- 
gion near the Fermi surface. It is natural then to assume that in /(k, k') the quasiparticle 
momenta lie on the Fermi surface. The function /(k, k') therefore only depends on the 
relative angle d. 

It is convenient to introduce the Fermi velocity and the effective quasiparticle mass 

by 

de(k) I * kp 



and to write 



/(k,k') = ^) (2.4) 

7T Z K 



where F($) is a dimensionless function which can be expanded in terms of Legendre poly- 
nomials, 

= X!( 2Z + 1 ) Fl P i( cos #) ( 2 - 5 ) 
i 

In the right hand side of (2.4) the prefactor is equal to the density of states on the Fermi 
surface. According to [18], the quasiparticle mass is related to the bare mass via 

^ = 1 + ^ (2.6) 
m 3 

An interesting property of the Fermi liquid at zero temperature is its compressibility, 

dP_N^_ 

U ~ dp~ mdN [ZJ) 
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where P is the pressure, p is the density and N is the number of particles. The value of u 
also defines the speed of (normal) sound in the Fermi liquid. According to [18], it is given 
by 

u 2 = (1 + F ) = V 4 (1 + F ) (l + ?f) (2.8) 



3mm* 3 \ 3 

where in the second equality we used (2.3) and (2.6). Another interesting property of the 
Fermi liquid is the existence of a massless excitation in the limit of zero temperature. This 
excitation, called zero sound, corresponds to the deformation of the shape of Fermi surface. 
The speed of zero sound uq is bounded from below by, and is generally proportional to, the 
Fermi velocity vf- To compute the ratio s = uq/vf one is instructed to solve the following 
integral equation 

(s - cos (p) = cos / F(tf>(tf', (ff) — (2.9) 

J 27T 

In eq. (2.9) the integral is over the solid angle and function ip) parametrizes the 
displacement of the spherical Fermi surface. 

In particular, when F($) contains only the zeroth Fo and first F\ harmonic [19] one 

has 

£ log fi±i)_ 1 = + 1 - (2 . 10 ) 

where s is the ratio between the speed of zero sound and the Fermi velocity 2 . In the limit 
of noninter acting quasiparticles, the speed of zero sound is equal to the Fermi velocity. In 
other words, s = 1 and Uq/u 2 = 3. [This result is independent of the particular form of 
F (•&).] Another interesting case is the limit of strongly interacting quasiparticles with Fo 
and F\ considerably larger than unity. From (2.10) we infer s 2 ~ F ° g Fl which combined 
with (2.8) leads to 

ir 

For generic and large Fi one can argue that the ratio between the speed of zero and first 
sound goes to a constant in the limit of strong coupling. This is because for a generic 
function F($) ~ F the right-hand side of (2.9) scales like F and therefore s ~ F. But 
according to (2.8) u ~ F as well, which justifies our assertion. 

This argument may fail for some non-generic values of Fi. In particular, when only 
Fo is non-vanishing and large, eq. (2.10) implies s ~ a/Fi/3. But in this case (2.8) again 
leads to uq/u ~ 1. 



2 Here zero sound velocity refers to the velocity of the mode for which v{-d, ip) is isotropic in 
the plane perpendicular to its momentum. 
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3. D3-D7 System 

In this section we review the D3-D7 system at zero temperature with nonvanishing 
chemical potential. Consider N c D3 branes extended along x , . . . , x 3 directions and Nf D7 
branes extended along x°, . . . , x 7 directions. Separating these branes in the x 8 — x 9 plane 
gives a mass to the hypermultiplet whose degrees of freedom are massless excitations of 
the fundamental string stretched between the D3 and D7 branes. Taking the near horizon 
limit, N c ^oo, g s ^0 with A fixed but large and Nf <C N c we obtain the AdS$ x S 5 geom- 
etry with N c units of RR five-form flux and D7-branes propagating in this background. 
Studying the holographic dual of the D3-D7 system with nonvanishing chemical potential 
involves analyzing the DBI action for the D7 branes embedded in AdS$ x S 5 with the 
electric field flux on their worldvolume turned on. The value of the chemical potential is 
equal to the asymptotic value of the gauge field on the brane. 

In the following we review the results of [16] where analytic results for the zero 
temperature limit of the D3-D7 system were first obtained. 3 The background metric can 
be written as 

ds 2 = (-dt 2 + dx l dxi) + (-0 2 (dr 2 + r 2 dn 2 + dR 2 + R 2 d(j) 2 ) (3.1) 

where t is time, i = 1, 2, 3 denote the spatial directions along the D3-brane, and dQ 2 is the 
metric of a unit three-sphere within S 5 . In eq. (3.1) p is the radial coordinate transverse 
to the D3 branes which is further expressed as 

p 2 = r 2 + R 2 (3.2) 

In this parametrization R and are the polar coordinates in the x 8 — x 9 plane, and U(l) 
symmetry allows one to set = const. Hence, the brane embedding is specified by a single 
function R(r). The induced metric is then 

ds 2 D7 = (|) 2 (-dt 2 + dx'dxi) + (-0 2 [(1 + R'{r) 2 ) dr 2 + r 2 dn 2 } (3.3) 

where R'(r) = The DBI action for this configuration is now relatively simple. Warp 
factors drop out and we have 

S D7 = -M J drr 3 ^Jl + R'(r) 2 -A' (r) 2 (3.4) 
3 Note that our notations are slightly different from those of [16]. 
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where Af = 



1 1 



N c Nf. Note that we divided by the volume of R 1 ' 3 therefore (3.4) is 



(2tt) 4 2A j 

actually an action density. Moreover, we fixed the gauge by choosing A r = and rescaled 
A according to: 

A ^27iA (3.5) 

Here and in the rest of the paper we set the string length to one. Given that (3.4) contains 
only first derivatives of the fields R(r) and A (r) there are two conserved charges 



R' 



A' 



= —c 



^l + R' 2 -(A' r ' ^l + R' 2 -{A' Y 

One can re-express R'(r), A' (r) in terms of c and d 



= d 



(3.6) 



R' = 



A' = 



Vr 6 + d 2 -c 2 ° Vr 6 + d 2 - c 2 



(3.7) 



The bare quark mass m and chemical potential \i are related to the asymptotic values of 
R(r— >oo) and Ao(r-^oo) as 



R(r^oo) = m = 2ivm, A (r^oo) = [i = 2nfi 



(3.8) 



Possible phases of this system are classified by the values of integration constants c 
and d. For c = d = 0, both R(r) and A (r) are constant. In this phase the condensate and 
the charge density vanish and the values of of R and A correspond to the quark mass and 
the chemical potential respectively. This is the zero temperature limit of the " Minkowski" 
embedding. Solutions of (3.7) for all other values of c and d satisfying d 2 — c 2 > are 
the zero temperature limits of the finite temperature black hole embeddings. They are 
characterized by the boundary conditions 



R(r = 0)=0 A (r = 0) = 



(3.9) 



and nonvanishing values for both the condensate and the charge density. Equation (3.7) in 
this case yields 



1 



R(r) = -c {d 2 - c 2 ) 



2\"l/3 



1 



A (r) = -d {d 2 - c 2 ) 



2\-V3 



B 



B 



r b 11 

r 6 + rf 2 -c 2; 6'3 
r 6 11 



r 6 +d 2 -c 2 ' 6' 3 



(3.10) 



where B r6+ ^_ c a ; g , 3 denotes the incomplete beta function while the values of c and d 
are related to the physical variables m = ^- and fj, = #- through 






Eq. (3.11) can be inverted to yield 



m = C7 -^ 2 -c 2 )-^ Ji = d 1 - 1 *(d 2 -c 2 )- 1 * k° F = fl 2 -m 2 = 7 -§(rf 2 -c 2 )^ (3.12) 



Combined with d 2 — c 2 > 0, (3.12) shows that black hole embeddings are only realized 
when the chemical potential is greater than the bare quark mass, [i > m. 

The phase structure of the D3-D7 system at T = has been analyzed in [16] with the 
following results. When // < m there is only one configuration available, the Minkowski 
embedding. In this phase, the system is characterized by zero condensate and charge 
density. For [i > m the configuration with brane falling to R = r = 0, the black hole 
embedding, is thermodynamically preferred. The two phases are separated by a quantum 
phase transition at ji = m. The black hole embeddings are distinguished by the non- 
vanishing values of both the condensate and the charge density. 



where c and d are related to m and [i via (3.11). Note that the expectation values in 
eq.(3.13) and, in particular, the charge density, vanish at the point of phase transition 
[i = m. In the rest of the paper we will be concerned with the black hole embeddings. 

4. Zero sound 

We will be interested in finding the massless excitation in the D3-D7 system at strong 
coupling. This requires analyzing linearized equations of motion which follow from the 
DBI action for the D7 brane propagating in AdS§ x S 5 . In this Section we only quote a 
few key results. For technical details the reader is encouraged to consult Appendix A. 

The DBI action can be written as 4 



The full action contains a Chern-Simons term. It can be shown that it does not contribute 
to quadratic order in the fluctuations for the field configurations considered below. 



(O m ) = 2ncM (J ) = 2ndAf 



(3.13) 




(4.1) 
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where Q = G + T is the sum of the induced metric and gauge field strength. Here we 
will only be interested in perturbations which are independent of the coordinates on S 3 . 
We therefore expand the D7 brane embedding coordinate R = R + £, the gauge field 
Aq = A^ + A , Ai = Ai and the gauge field strength jF 0r = + F 0r , J 7 ^ = Fij, where 
i,j = 1,2,3 and the superscript (0) denotes the background values of the corresponding 
fields (whose profile follows from the results of the previous section.) In addition, we set 
Aq — 0,V 6 G S 3 and choose the A r = gauge. To expand the DBI action to quadratic 
order in fluctuations one can make use of the formula 

Vtetg = VdetGo + ^ (trg-^g) 2 - ^{g^sgg^sg^j (4.2) 

where we write g = go + dg and 5g stands for the fluctuations. In (4.2) there is of course a 
linear term, but it vanishes by equations of motion. It is convenient to keep the fluctuations 
of the brane embedding £ (but not derivatives of £) as a part of (/o- Both the leading term 
and a term linear in fluctuations do not contain £, which means it does not appear at 
quadratic order in fluctuations. 

The action for fluctuations is given by [see also eq. (A.l)] 



S D7Jl = -\M J d 4 x J drg(r) 



i<3 



+h(r)(d r 2 - — (dot) 2 + L 4 / 4 (r) ~ 2L% £ F 0i (d^) - 2/ 7 (r)F 0r (<9 r £) 

(4.3) 

where the functions i = 1, . . . , 7 can be found in (A. 2) and (A. 3). It is convenient to 
use the momentum space representation, 

A M (x",r) = J 0^e ik ^A M (k»,r) (4.4) 



and 



d 4 k ~ 

(U".r)= / e **M*^(fcM, r ) (4.5) 



We can also choose the direction of wave propagation by setting k^ = (— uj, 0, 0, q). It is 
important to work with the gauge invariant combination, 



E = qA + uA\\ 
8 



(4.6) 



We can use Gauss' law (equation of motion for A r ), 

qA\+uhA! -uf 7 l' = (4.7) 

to express A' (or A'^) in terms of E' and Here prime denotes differentiation with 
respect to r. In these equations and in the following we define q = ql? = q\f2\ and 
uj = uL 2 = uV2X. Equations of motion are in general quite complicated. In particular, the 
equations for the longitudinal components of the gauge field are coupled to the equations 
for the embedding fluctuation £. It turns out, however, that they can be solved analytically 
in the two overlapping regimes, much like in [14]. 

As explained in Appendix A, in the vicinity of the horizon equations for E and £ 
decouple: 

E+ -E + n 2 E = (4.8) 

z 

't+-J+n 2 z = o (4.9) 

where O 2 = ui 2 d ~^ and the dot denotes differentiation with respect to z = 1/r. The 
solutions of (4.8) and (4.9) are 

E(z) = A , m = B (4.10) 

z z 

The positive exponent is singled out by incoming boundary conditions at the horizon [20- 
23], [14]. In the limit of small u, Qz < 1 (4.10) becomes 

E(z) = itlA + -, £ = itlB + — (4.11) 

z z 

On the other hand, for sufficiently small ui and q equations of motion simplify again: 

z hf 5 q 2 - u 2 f 1 ) \ hf 5 q 2 - u) 2 f x ) 

(4.12) 

z hj5 T - u 2 fi 

where Z = g(z)E(z), Y = qg(z)£(z), and h{z) = (1 + (d 2 - c 2 )z 6 )/(l - c 2 z 6 ). The solution 
of (4.12) is given by 

E = 3- [d[(l - c 2 z % )^ - (1 + (d 2 - c 2 )z 6 )u 2 } + ^cdz^q 2 ] 

(1 + (d 2 - c 2 )z 6 )* ' ' x 

x. ( 4 - 13 ) 

£ = q - 3- [dcdz 6 - C 2 (l + d 2 z 6 )] 

(l + (d 2 -c 2 )z G )z 



where C\ and C 2 are arbitrary integration constants. Eqs. (4.13) can also be obtained 
directly from eq. (A. 6) by neglecting all non-derivative terms, performing trivial integra- 
tion, and using Gauss' law to express A' in terms of E' and Eqs. (4.13) can in turn be 
integrated to yield 

E = C + f — r [Ci[(l - c 2 x & )f - (1 + (d 2 - c 2 )x 6 )u 2 ] + C 2 cdx 6 q 2 } 

Jo (1 + (d 2 - C 2 )x^ " 

l=Co+ T ^ 3- [C lC dx 6 - C 2 (l + d 2 x 6 )} 

Jo (i + (dp - c 2 )x 6 y 

(4.14) 

The integrals can be expressed in terms of hypergeometric functions. We do not need this 
representation, since only small and large z asymptotics are important for our purposes. 
In particular, near the boundary 

E = C + O(z 2 ), t = Co + 0(z 2 ) (4.15) 

This implies that the spectrum of quasinormal modes can be obtained by requiring Co = 

c = o. 

In the region z 3> 1 the leading asymptotic behavior of (4.14) is a constant whose 
value we can infer by performing integration from z = to z = oo. The subleading 1/z 
term can be extracted from the expression for the derivatives (4.13). The results are 

E = Co + ftxd + b 2 C 2 + ^ + ^ (4.16) 

z z 

where 

c 2 ^ + (d 2 - c 2 )u 2 cdq 2 

a\ = o , a 2 = s- (4.17 

{d 2 -c 2 )i (d 2 -c 2 )l V ' 

and 

h = -^^-[(3c 2 -d 2 )f + 3(d 2 -c 2 )u 2 }, b 2 = -2cdf ^Pl}P (4.18) 
Likewise, 

£ = Co + 6^ + b 2 C 2 + + (4.19) 

where 

~ cd „ d 2 , , N 

ai = 3-, a 2 = 3- (4.20) 

(rf 2 -c 2 )i (rf 2 -c 2 )i 
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and 



~ 2cdv a r ~ r r , , . 

6 = Ve; 4 V3; , b 2 = ^ — (3d 2 - c 2 4.21) 



We will focus on the quasinormal mode with the linear dispersion relation in the regime of 
small uj and q. In this case we can match the near horizon solutions (4.11) to (4.16) and 
(4.19). Requiring Co = Co = and neglecting terms of higher order in u and q we obtain 

hd + b 2 C 2 = 0, b 1 C 1 + b 2 C 2 = (4.22) 

These equations can be simultaneously solved with nonvanishing C\ and C 2 whenever 
bib 2 — b 2 bi = 0. Using (4.18) and (4.21) we arrive at 

d 2 - r 2 

~2 2-2 2 u c (a oo^ 

u = u q , u = M2 _ ^ (4.23) 
Using (3.11) we can write the expression for the speed of zero sound as 

2 2 

This is the main result of this paper. In the relativistic limit |j»mwe recover the result 
of [14], Uq = 1/3. In the vicinity of the phase transition, where Sfi = fi — m m, the 
speed of zero sound vanishes as Uq ~ 8ji/m. As discussed below, precisely such a behavior 
is expected from the Landau's theory of Fermi liquids at strong coupling. 



5. Discussion 

Before discussing the significance of (4.24) let us make a comment regarding the 
identification of the gap less mode with zero sound. At first sight, our calculations imply 
that the dispersion relation is modified by a term of the type — iq 2 . In particular, in the 
simpler case of vanishing mass considered in [14], the dispersion relation looks like 

w = ±--i-r (5.1) 
<j as 

Note that in [14] the —iq 2 behavior of imaginary part of the pole was important for iden- 
tification of the massless mode with zero sound. In the derivation of dispersion relation it 
has been assumed that uz 1 and rf3z»l regions overlap. More precisely, the solutions 
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(4.8), (4.9) and (4.13) are valid up to terms 0(d- 2 z~ 6 ) and 0(w 2 z 2 ) respectively 5 . This 
seems to imply that (5.1) is only corrected at higher order in q. 

A number of quantities in the Fermi liquid setup can be computed if the particle 
density is known as a function of Sfi = \x — m. In principle, we can compute the charge 
density using 

It would be equal to the quark density n, if the quarks were the only degrees of freedom. 
This would allow one to determine fcp as a function of \i. However the situation is more 
complicated, as both the fermions and the scalars in the fundamental hypermultiplet are 
charged under U(1)b- At vanishing coupling the system is unstable due to the condensa- 
tion of bosons. However the lagrangian contains a term which is quartic in scalar fields. 
Hence, at large A the condensation of bosons can be stabilized, and indeed the solution for 
the brane embedding does not exhibit instabilities. At large coupling one can use the holo- 
graphic dictionary to compute the expectation value of the operator which contains the 
scalar bilinear. And, as noted in [16], the dependence of this condensate on 5[i is consistent 
with the quartic potential for the scalars and a quadratic term of the form (m 2 — n 2 )qq* ■ 
In summary, using (5.2) to determine kp would yield hp ~ (<i/A)i ~ (m 2 5/u/A)i but 
this result is unreliable due to contribution of the boson condensate to the charge density. 
However the expression for p can still serve as an upper bound on n. 

We can deduce the critical behavior of various quantities using the results reviewed 
in Section 2. Writing (2.7) in the form 

« 2 = ^ (5.3) 
m on 

and noting that n vanishes near the critical point faster than dfi, we infer u 2 ~ 5p> and, 
consecutively, vf ~ 5(1^ . This implies a dispersion relation of the form e(k) = k 2 /2m* + 
0(k 3 ). Hence, in the vicinity of the critical point the density of quarks behaves like 
n ~ k F ~ 5fi%. This, in turn, leads to 

u 2 ~ ™± (5.4) 
3m v ; 



We thank Andreas Karch for very useful discussions on this issue. 
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As reviewed in Section 2, the speed of zero sound is proportional to Fermi velocity as 
well, and must therefore vanish at the critical point as u ~ This is precisely what 

eq. (4.24) implies in the limit Sfi <C m: 



« ^ (5.5) 



m 



Hence, we observe that our string theoretic computation reproduces (in a rather nontrivial 
fashion) the critical exponent predicted by the phenomenological theory. 

What about the coefficient in (5.5)? As reviewed in Section 2, in the limit of 
noninteracting fermions = vp or, in other words, (u^) 2 = 25^/m. This implies 

/u^ = \/3. However this ratio is expected to be modified when the coupling is strong. 
Indeed, our result uq/u = a/3/2 differs from the value at weak coupling by a factor of 
This is not surprising since, as discussed in Section 2, the ratio of the velocities of first 
and zero sound goes to a constant which is 0(1) for generic strong interaction between 
quasiparticles. 
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Appendix A. Fluctuations 

The action is 



S D7Jl = - l -M J d A x J drg(r) 



E F ir~ f^ F 0r ~ L 4 h(r) £ Fl + L 4 / 3 (r) £ F* + 



i<3 



L 



4 



+h(r)(d r 2 - —(doO 2 + L 4 h(r) 5>,£) 2 - 2L 4 / 6 (r) £ F 0i {d^) - 2f 7 (r)F Qr {d r i) 
Po 



(A.l) 
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where the functions f t with i = 1, • • - 7 are defined in terms of the background fields as 
follows 

- 3 1 + R' 2 , , , 1 + R' 2 



g(r) 

Mr) 



1 + R' 2 - F^ 2 



Mr) = 



1 + R' 2 - F n (0)2 



Or 



Mr. 



1-F, 



(0)2 



Or 



1 _ P(°) 2 

Mr)= ° r 



1 + R' 2 - F ( r 0)2 



(o) 



Mr) 



Pt 



Mr) 



R'r F (0) 



1 + R' 2 -F^ 2 



Using (3.7) we can rewrite these functions as 



g (r) = Vr 6 + d 2 -c 2 /i(r) = l + 



d 2 



Mr) 



Mr) 



r 6 + d 2 — c 2 Pq 



Mr) 



Mr) 



cd 



^ + d 2 - c 2 p A 



r b -(? 1 
r 6 + d 2 — c 2 p^ 

Mr) = '4 



r b + d 2 1 
r 6 +d 2 -c 2 tf> 

Mr) = l-~ f . 



(A.2) 



(A.3) 



with po(r) given by (3.2). It is useful to note that not all of these functions are independent 
from one another. For instance 



Ji — fzfi fe — ^3/7 — fzfs 



(A.4) 



These identities will be particularly useful in expressing the equations of motion in a 
compact manner. The field equations are more conveniently expressed in the momentum 
space representation 



and 



A M {x»,r) 



£(x",r) 



d^k 
d 4 k 



e^ x ^(^,r) 



A M {k^r) 



(A.5) 
(A.6) 



Further choosing k^ = (—u>, 0, 0, q) we find 

d r [gMd r A ) - gMdrO] - mgh\ - ^9/2^0 + q 2 gkt = 

d r g{d r A\\) + ojqgf 2 A + u) 2 gf 2 A\\ - toqgf 6 ^= 
(),. , 

d r [g(d r A ± )] + g[u 2 f 2 - ?f 3 ]A ± = 



9MdrO + gMdrAo) + u 2 ^- fgM - qZgfcA {l - fghAo = 

Po 



(A.7) 



14 



Here we have absorbed the 't Hooft coupling constant A into the variables uj,q defined as 
q = qL 2 = q\/2X and uj = ujL 2 = ujV2X. To fix the residual gauge invariance we impose 
Gauss law 

qA\+uhA' Q -uf 7 l' = (A.8) 
Choosing the gauge invariant combination 



E = qA + uA\ 



(A.9) 



and using (A.8) to express the first derivative of Aq in terms of E' and £' 

',2 



Ar 



q 



q 2 - Z 2 h 



UJ 



E' - — M' 
Q 



leads to 



t:" + [ 9 - + f i + yt\„ ) E> + / 3 (u/ 2 /i - f)E- 



9 fi T ~ w f 



fl 

h q 



9' , fh , u 2 f'i 



X" + [— + -j- 

\9 h 

7/2 



7a 



u 2 h 



X' + fs(u 2 h - (f)X 



= 



feq 2 - hu 2 



9' , # , Z 2 f[ 



E n , o_ + <j_ 



9 fl Q 2 ~ u 2 fi 



E' + fs^h-q 2 ) 



E+ 



+ X" + 



9' . ft , Z 2 f{ 



UJ 



2 h! \ 



where 



g f s q 2 — uj 2 fi q 2 -uj 2 h) 

1 Pi + hh r 6 + (d 2 - c 2 ) 



X' + f 3 (Z 2 f 1 -q 2 )X = 



h 



r 6 — c 2 



(A.10) 



(A.ll) 



(A.12) 



and we defined X = q^. Performing a change of variable from r to z = £ (A.ll) is recast 
into 



E+\ 2 - + 9 - + ^ + 



UJ 2 h \ i, , h,~2 



E+^h-^E 



z 9 fi Q 2 — u 2 fi 



fl 
fi 



z 9 fl Q ~ w fi 



X + ^h-^X 



fi Q 2 _ 

f 5 q 2 - hu 2 



E+[ 2 - + 9 - + t + 



u 2 h 



z 9 fi T ~ u 2 f 



(A.13) 



+ 



2 9 , h , ^ 2 /i 



x + 4(5 2 A-^)x = o 



2 # /5 Q 2 - ^ 2 /l ? - ^ 2 /l 
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where dots indicated differentiation with respect to the variable z. Observe that (A. 13) 
reduces to two independent equations for E and X near the boundary since the prefactors 



h 



cdz^ 



fr Q 2 



cdz t 



s 2 



h l + d 2 z 6 hqt-hu 2 l-^z^-u 2 
vanish when z — > 0. We are then left with 



(A.14) 



E(z) — -E(z) + (uj 2 — q 2 ) E(z) = X(z) - -X(z) + (u 2 - q 2 ) X(z) = (A.15) 

2/ £ 



The general solution of this equation can be written in terms of Bessel functions of the 
first kind as 

E(z) = AM [z^w 2 - q 2 ] + &\J\ \zV~u 



' uj 2 — q 2 



X(z) = A 2 Y 1 
Expanding around z = yields 



<^u 2 - Q 2 ] + B2J1 [zV^ 2 ~ Q 2 ] 



(A.16) 



E(z) ~ Ax 
X{z) ~ A 2 



l + \{u 2 -f) (l-27-ln 
l + l^-q 2 ) (l-27-ln 



^-q 2 2 
- z 



uj 2 — q 2 2 
: Z 



+ B lZ 2 
+ B 2 z 2 



(A.17) 



This identifies the constants Ai and Bi for i = 1,2 with the coefficients of the non- 
normalizable and the normalizable mode respectively. 

In the vicinity of the horizon on the other hand notice that 



9 3 1 h 6 

r^j r^j 

g z (d 2 — c 2 )z 6 fi z 
h^Q fr = 6 u 2 h 
h z fr z 



Z 2 h 



6 



q 2 - uJ 2 fi 



6 d 2 



uj 



q 2 -u 2 h zc*z 6 q2 + d2 - c2 Z 2 



(A.18) 



whereas 



(u 2 h - g 2 ) ^ 



d 2 



(A.19) 



since Rq(z) behaves for large z as ^ d Z_ c2 \ • It follows that eq. (A. 13) reduces to 



E +-E + Q 2 E - - 1 

z d 



x + -x + n 2 x 

z 



Q 2 



x + -x + n 2 x + -- ., , 

z cq 2 + d 2 C uj 2 



E + -E + Q 2 E 

z 



(A.20) 
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where we defined O = cuW ^ c ■ Multiplying either of the equations with appropriate 
factors and adding them up yields the following system of decoupled equations near the 
'horizon' 

e + -E + n 2 E = o x + -x + n 2 x = o (A.21) 

with solutions 

piYLz p iYLz 

E (z) = A- — X(z) = B- — (A.22) 

z z 

Note that the incoming wave boundary condition singled out the positive exponent in 
(A.22). 

Eq. (A. 13) can be simplified even further. 

E+ {-z + - g + W 5 ¥^) E ^ { " fl - q)E -{hh¥^h) x = 
x + ( 2 - + i + A : ?z«LU + A (i3S/i _ ?)x+ |'A = ?.U = 

\2 r - y ^ 4 \hf5q 2 - uj 2 h J 

(A.23) 

To arrive at (A.23) multiply either of the equations with the relevant factor and add it to 
the other one. Note that the equations for E and X are now coupled only through first 
derivative terms of the fields. This form is particularly useful when considering the region 
VLz <C 1. In this case, terms without derivatives of the fields E and X can be neglected to 
yield: 



2 , 1 Q 2 fi 




z + hf 5 q 2 -u 2 f x 

(A.24) 

2 hf->q-^ 2 Ji) yhf 5 q 2 -u 2 f 1 

where Z = g(z)E(z) and Y = g(z)X(z). Solutions for Z and Y are given by 

_ fCi [g 2 -u 2 {l + (d 2 - c 2 )z 6 )} + q 2 ^ [-q 2 (1 - c 2 ^ 6 ) + u 2 (l + (rf 2 - c 2 )z 6 )] 
Z ~ cdq^z 2 [1 + (d 2 - c 2 )z Q ] 



z 2 1 + (d 2 - c 2 )z G l + (d 2 -c 2 )z 6 

(A.25) 

where the integration constants C\ and C 2 are related to the ones appearing in the main 
text through 

C x = -C20 2 C 2 = -d (dC 2 - cCi) q 2 (A.26) 

Eq. (A.25) can be easily integrated to yield an expression for E and X in terms of 
hypergeometric functions. 
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Appendix B. The matching technique in detail. 

Here we investigate in detail the regions of applicability of the matching technique 
used in section 5. Let us focus on understanding the precise conditions under which eq. 
(A.23) reduces to either of eqs. (A.21) or (A. 24). 

Observe that in the vicinity of the horizon, rather when z is large enough so that 
z ^> (d 2 — c 2 )~^ > the following simplifications occur 

g 3 1 1 cfh 6 d 2 1 



g z (d 2 — c 2 )z 6 hfs q 2 — iv 2 fi z d 2 — c 2 ]_ _ ^ ^ljo 



7f2 ~2 c J2 1 — 



fs q 2 -uj 2 6 c 



hhq*-&h~ zd 2 -c 2 1 _Zl d 2 z e ( R1 ) 



g 2 



f 7 q 2 -u 2 6 cd 1 ~^2 f 7 q 2 6 cd 1 



V5?-w 2 /i z d 2 - c 2 i _ ^ld 2 z 6 hf b q 2 -uj 2 f 1 z d 2 - c 2 i _ ^d 2 z & 



q 2 q 2 



Moreover, given that R(z) for large z behaves like R(z) ~ ^Ji-c 2 z 

h(~2 t ~2x ^(^-C 2 ) 1 (u 2 2 6 \ 

{u h-q)- to. z - 1 ( B - 2 ) 



z* y J1 ^ ' ~ d 2 d 2 z 6 yd 2 
and if we additionally assume that z 3> d~^ (^i^j 3 eq. (A.23) reduces to (A.21). 

E + -E + n 2 E = X + -X + Q 2 X = (B.3) 

z z 

with O defined to be O 2 = uj 2 d ~ 2 C . The solution of (B.3) is given by 

iClz iViz 

E(z) ~ X(z) ~ (B.4) 



When furthermore z <C q the solution becomes 

E(z)=Ain+- E(z) = Bifl+— (B.5) 

z z 

with A, B overall multiplicative constants. In summary, (B.5) consistently describes 
the behavior of the solution of eq. (A.23) as long as z lies within the annulus 
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Max 



(d 2 -C^-^d-s 



<C z <C o- This in turn implies that we are effectively 



exploring the region of parameter space in which 



Max 



(cf -c 2 )-t,cH (I 



1 



(B.6) 



Now let us investigate the behavior of (A. 23) in the regime Qz <C 1 first. Note that 
terms without derivatives always appear with the prefactor |§ (p 2 fi — q 2 ) as compared 
to the two derivative terms. Moreover, this prefactor is given by the difference of two 
monotonic functions. Its absolute value is therefore bounded above by the absolute value 
of their sum. The latter is a monotonic function which tends to a constant for both large 
and small z. In particular, 



h(z) 



( w h + 9 ) 



l~2 i ~2| 



z ^> Max 



z < Mm 



(rf 2 -c 2 )- 



(B.7) 

It follows that terms without derivatives can be neglected as long as they are small com- 



pared to the two derivative terms. This means that for z <C Min 
(A. 23) consistently reduces to 



fi' Z 



1 + 



eq. 



Z + 



2 1 

+ 



eh 



z hfaq 2 - uJ 2 fi 
z hfbQ 2 - u 2 fx 



Z - 



Y + 



fr q z -Z 2 
hfc q 2 - u) 2 f\ 

fr ^ 
hh q 2 - u 2 fi 



Y = 



Z = 



(B.8) 



with Z = g(z)E(z) and Y = g(z)X(z). 

As previously explained, we would like (B.3) and (A. 24) to be simultaneously valid. 
This implies that we can investigate (A. 13) in the following region of parameter space 



Max 



{d 2 -C 2 )-^d- L z 



UJ 



< Min 



1 1 



1 + 



Q 2 



(B.9) 



It is manifest from (B.9) that the analysis of section 5 covers the region of small uj, q. 
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